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ABSTRACT 

Electromagnetic processes associated with a charged particle moving in a strong circular 
magnetic field are considered in cylindrical coordinates. We investigate the relation between 
the vacuum curvature emission and Cherenkov emission and argue that, for the super luminal 
motion of a particle in the inhomogeneous magnetic field in a dielectric, the combined effects 
of magnetic field inhomogeneity and the presence of a medium give rise to the synergetic 
Cherenkov-curvature emission process. We find the conditions when the operator relations 
between electric field and electric displacement in cylindrical coordinates may be approximated 
by algebraic relations. For nonresonant electromagnetic wave, the interaction with particles 
streaming along the curved magnetic field may be described in the WKB approximation. For 
resonant waves, interacting with superluminal particles we use a plane wave approximation to 
compute the local dielectric tensor of a plasma in a weakly inhomogeneous magnetic field. We 
find in this approximation the polarization of normal modes in the plasma, Cherenkov-curvature 
and Cherenkov-drift emissivities and growth rates. 



Subject headings: stars:pulsars-plasmas-waves-radiative transfer 



Introduction 



Despite much theoretical effort, there is still no widely accepted explanation of how pulsars emit the 
high brightness radio emission by which they were first discovered over thirty years ago. Difficulties with 
"antenna" mechanisms have led to a resurgence of interest in "maser" mechanisms where an inverted 



population of electrons and/or positrons amplifies an outgoing wave mode (e.g. Melrose 1995). Two such 
maser process that are particular promising are the Cherenkov-drift mechanism (which is Cherenkov- 



curvature emission with a drift ( Lyutikov, Machabeli fc Blandford 1998 )) and cyclotron-Cherenkov emission 
at the anomalous Doppler resonance ( Kazbegi et al. 1991 ). This paper is concerned with elucidating the 
physics of the Cherenkov-curvature and Cherenkov-drift processes and developing a new mathematical 
description of it using cylindrical coordinates. Using this approach, its close relationship to the pure 
Cherenkov emission processes becomes clear. 

Pulsar radio emission is believed to originate on the open magnetic field lines that trace a path from 
the neutron star surface to the light cylinder at r = c/fl and beyond. The field geometry is complex 
close to the star at r « i?* ~ 10 km and near the light cylinder. However, it should be primarily dipolar 
for i?„ <C r <C c/fi. The characteristic length scale, essentially the radius of curvature of the field lines, 
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Rb ~ (cr/fi) 1 / 2 , is very large compared with the wavelength of interest, A, and so we should be able to 
adopt a WKB approach following following individual wave packets as they propagate out through the 
magnetosphere. However, in order to compute the local emission and absorption it is necessary that the 
magnetic field remains curved. A simple model problem, that retains the essential ingredients, comprises a 
set of circular concentric cylindrical magnetic surfaces. (Fig. |l|). We show, that it is possible to ignore the 
radial variations in the strength of B in computing the local interaction, provided that this variation is not 
very rapid. In this model problem the plasma circulates continually along the circular trajectories. 

It is generally presumed, that electrons and positrons are in their ground gyrational state and follow 
the curved field lines. This is a good approximation for the bulk of the plasma, which we will suppose 
travels with a Lorentz factor 7 P »s 100. However, we believe that there is also a population of ultra energetic 
particles with 7 = 10 5-7 and, in the outer magnetosphere, these will experience a curvature drift relative to 
the bulk plasma. This has the interesting consequence that their main electrodynamic interaction is with 
the waves propagating in the bulk plasma at a finite angle to the magnetic field B. It is these waves, that, 
we assert, become the high brightness electromagnetic waves that escape from the magnetosphere. 

When considering the wave amplification in this problem, we assume that fluctuating currents, present 
on the field lines at radii smaller than the radius of the considered region, produce electromagnetic 
fluctuations that are subsequently amplified due to the interaction with resonant particles (Fig. |Tj). This is 
a convective type instability, when a wave is amplified as it propagates through an active medium. 



A first attack on this problem was made by Blandford 1975 In this paper the energy transfer from a 



plane infinite electromagnetic wave to a single electron moving with ultrarelativistic speed along a curved 
trajectory was found to be always positive, independent of the distribution function and so there is no 
possibility of wave growth. The subsequent works ( Zhclcznyakov fc Shaposhnikov 1979 , Lou fc Melrose 19*9^ , 



Melrose 1978j ) basically followed this approach, which emphasizes the analogy between curvature emission 



and conventional cyclotron emission. This approach, though formally correct, has limited applicability 
and ignores two important features of the emission mechanism. The first is that, in adopting a plane 
wave formalism, the interaction length for an individual electron, sa Rs/jb, was essentially coextensive 
with the region over which the waves could interact with any electron. This approach precludes a strong 
amplification under all circumstances because the wave would have to grow substantially during a single 
interaction. The second problem was that a dispersion of the waves was neglected. We address the first 
shortcoming by expanding the electromagnetic field in cylindrical waves centered on r — 0, and the second 
explicitly by considering general plasma modes. 



In a separate approach developed by Bcskin, Gurevich fc Istomin 1986 an attempt was made to 



incorporate the collective effects of a plasma in an inhomogencous magnetic field. The fundamentals of that 



approach have been seriously criticized (Nambu 1989, Machabeli, 199?;) 



One of the key new elements in our approach is that we expand the electromagnetic fields in cylindrical 
coordinates and consider resonant interaction between a particle and these modes. A common procedure 
in calculating the energy emitted by a given current distribution is to find the power emitted into a given 
normal mode of the medium and sum over all the modes. A power emitted into a normal mode by a given 
charge distribution is proportional to the square of the expansion of this current in normal modes. In a 
homogeneous or weakly inhomogeneous stratified medium considered in Cartesian coordinates the normal 
modes are plane waves (with slowly changing parameters), so that the power emitted into a given normal 
mode turns out to be proportional to the Fourier transform of the current. To find a similar expression 
for a single particle emissivity into a cylindrical mode in an inhomogeneous medium is difficult problem 
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for two reasons: there is a complicated, radius dependent, form of the vector cylindrical waves and the 
dielectric response in an inhomogeneous medium is nonlocal. We address these complications by finding the 
conditions, when the interaction of a particle with cylindrical waves can be approximated by an interaction 
of a particle with plane waves. We find two cases when this can be done. First, this may be a true 
approximation for thenonresonant modes. This is equivalent to the WKB approximation to the radial 
dependence of normal modes (mathematically, this corresponds to the tangent expansion of Bessel functions 
when argument is significantly larger that he order). In this case the response of a medium becomes local. 
Second, we find a particular case, when the resonant modes can be approximated as local plane waves. 
A resonant interaction of a relativistic particle with a cylindrical mode occurs near the point when the 
argument of Bessel functions is close to the order. The WKB approximation, or expansion in tangents, is 
not applicable in this case and we have to use the Airy function approximation to Bessel function, which 
has a plane wave approximation for the interaction of subluminous waves with the particles moving with 
speed larger than the speed of light in a medium This corresponds to the Airy function expansion argument 
being larger than the order. 

We should note here, that the electrodynamics of the interaction of a particle moving along the 
curved magnetic field with speed larger than the speed of light in the medium is quite unusual and can be 



considered as a new, Cherenkov-Curvature, emission mechanism ( Lyutikov, Machabeli & Blandford 199 



which differs from conventional Cherenkov, cyclotron or curvature emission and includes, to some extent, 



the features of each of these mechanisms. In an extension of ideas of Schwinger et al. 1976, conventional 
synchrotron emission and Cherenkov radiation may be regarded as respective limiting cases of \n — 1| — ► 1 
and B — » of a "synergetic" cyclotron- Cherenkov radiation, the Cherenkov-curvature radiation necessary 
includes effects of the magnetic field gradients. 

When the dielectric response of a medium to cylindrical wave is local it is possible to calculate a 
simplified dielectric tensor. Using this dielectric tensor we find the normal modes of strongly magnetized 
electron-positron plasma and show that a beam of particles propagating along the curved magnetic 
field can amplify the electromagnetic waves. The amplification occurs at the Cherenkov-drift resonance 
lu — k^v^ — k x Ud = (u> is the frequency of the wave, k$, k x and are the corresponding projection of 
the wave vector and velocity and Ud is a curvature drift velocity). The fact that this is a Cherenkov- type 
resonance immediately implies that the presence of a subluminous waves with the phase velocity smaller 
that the speed of light is essential. The presence of a drift provides a coupling between the the electric 
field of the electromagnetic waves particle's motion. Our estimates show, that this instability can grow fast 
enough to account for the observed pulsar radio emission. This instability may be regarded as new type of 
a curvature maser. An interesting and peculiar feature of this mechanism, is that, unlike with conventional 
curvature emission, the emitted waves have a polarization almost perpendicular to the osculating plane of 
the magnetic field. 

The overview of this work is the following. In Section ^| we analyze the properties of electromagnetic 
waves in cylindrical coordinates in vacuum and the interaction of the vacuum waves with a charged particle. 
In Section ^ we analyze in cylindrical coordinates the electromagnetic fields of a particle moving along a 
spiral trajectory using the dyadic Green's function for the vector wave equation, derive the emissivity of a 
particle in ground gyrational level into a cylindrical mode and rederive the curvature emissivity in cylindrical 
coordinates. This is followed in Section ^ by the generalization to a dispersive isotropic medium and the 
importance of waves with the phase speed less than the speed of light in vacuum is brought out. In Section 
^ we consider anisotropic plasma in infinitely strong magnetic field and derive the curvature emissivity 
using the Vlasov approach. In Section ^ we discuss the various regimes of the Airy function approximation 
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to the Cherenkov-curvature emission and find the conditions for the plane-wave approximation. In Section 
we investigate the features of the electromagnetic waves in the asymptotic regime z> 1 and in Section ^| 
we calculate the response tensor for a one dimensional plasma in a strong curved magnetic field taking into 
account the drift velocity. Finally, in Section || we analyze the polarization properties of electromagnetic 
waves in the plane-wave approximation and calculate the increment of the Cherenkov-drift instability. 



2. Vacuum Solutions 

2.1. Vacuum Normal Modes 

We expand the fundamental solutions of the wave equation 

u 2 

curl curl E -E = (1) 

in terms of Fourier amplitudes in x and <fi coordinates and time: 

f=oo „ I" dk 

E(r,t)= du —^E(r,k x ,v,u)exp{-i(wt-v^-k x x)}. (2) 

V — — oo 

The wave equation (Q) then takes the form 

iv d , „ . ., d „ k r v „ ( v 2 ( uj 2 



__ {rE ,) + ik x -E x - -^E, + [--[--kl\\E r = Q (3) 



d (\ d , „ A . d f E r \ k x v 



,2 



' " ^lE x ) + i ^l(rE r )-(4-t T )E l 



r dr \ dr J r dr 
Using the fact that in a wave the magnetic field is related to the electric field: 

curlE = — =i-B. (6) 

c at c 

we find the amplitudes E(r, k x , to), subject to condition of being finite at r = 

E^(r,k x ,u,w) = E^ f^Jl(Xr)e r - ^ J„(Ar)e^ + J,(Ar)e x ) (7) 

BC'lfr,^,^)^^) (^!l^Me r + zAj:(Ar)e^ (8) 

E^(r,k x ,u,uj) = E^ f^.J u (Xr)e r -J' v (\r)eA (9) 

B (t) (r,t,^) = f^(Ar)e r + ^J„(Ar)e - i- J„(Ar)«J) (10) 
V u Arw w / 



where A = y^uj 2 /c 2 — k 2 , J„(Ar) are Bessel functions, which satisfy the equation 

and e r , and e x are unit vectors along the corresponding axes, E^ lt > and EW are amplitudes of the modes. 
The superscripts denote polarization: E^*) is perpendicular to the coordinate surface e x and E^ lies on 
the coordinate surface e x . 



- 5 - 



2.1.1. Short Radial Wave Length Expansion of Vacuum Solutions 

For large orders v 3> 1, argument larger than the order Ar > v and not very close to it, Xr — v > v, we 
can use expansion in tangents of the Bessel functions: 



J^secC)~J tCOS KtanC-C) + T ), C > ^ 1/2 (12) 

ttv tan Q V 4/ 



In our case cos ( — v/ (Ar). 
Introducing 



kl = A 2 -, (13) 



the Eq. for the Bessel function may be viewed as a dispersion relation: 



K + kl + k 2 x = - (14) 

Then the expansion (|l^) reads 

J„(Ar) w \ -y— cos (±k r r- <f> v ) (15) 
V 7rfc r r 

where 0^ is insignificant phase shift. The limits of applicability of this expansion is that k r r 3> 1 and £ not 
very small. The normal modes in this limit are 

E< tt > (r, t) oc e -^-k x ,±k r r) (_ Mp^ _ *x** ^ + ^ 

E« (r, t) cx e-*^*-*-**^) f^e r - ^e/) (16) 



2J.i?. WTfB Solution (Short Radial Wave Length) 

We can obtain solutions (|l^) of Eqns (j^- 1|) using WKB procedure, when the effective radial wavelength 
is much shorter than the characteristic scale of the problem. Assuming that the solutions of Eqns (|^- 1|) 
have a form 

Eoc-L e * s , (17) 
yr or r 

we find 



S = ± J dry A 2 j — i (^V r — ^ arccos — ^ (18) 

with fc r defined by Eq. (|l4|). 

The dispersion equation, which is obtained by substituting Eq. ( p~7[ ) in Eqns ^|), takes a form 



where 



4?^ = (19) 



/ fc 2 + /c 2 fef&# + "2^" k x k r i^j- 

-k r kt-^ fc 2 +fc 2 + 3 ^-^ -k x k4, | I 20 1 

_l, jL I ifcx _fr fr. h 2 4- h 2 — 1 — u ' 
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For k r 3> 1/r we can drop the complex part of the tensor (which is not even Hermitian since k r is 



not a Killing vector). The dispersion equation det||e^|| = then gives again 



u 2 

k r + k x + — —j- (21) 

The Bessel equation ([ll]) may be viewed as a Schrodinger type equation for a wave functions of a 
particle moving in a two dimensional potential given by the v 2 jr 2 term in Eq. (|ll]). The characteristic scale 
of the Schrodinger-type equation ( |ll|) is r v = v/\. A point r = r v may be viewed as a classical reflection 
point. Solutions are exponentially decaying for smaller radii and have a wavelike structure at larger radii. 
The analogy between the Eq. ([ll] extends even further. When approaching the classical reflection point 
a common approach in quantum mechanics is to use the Airy function approximation to the solutions of 
the Schrodinger equation. This corresponds to the transition region in the Bessel function expansion, when 
relations (|l2|) are no longer valid and one should use Airy function expansion. We will discuss the various 
regimes of the Airy function expansion in Section ^. 

The WKB method breaks down when the radial wavelength becomes comparable to the inhomogeneity 
scale of the " potential" . For a given radius this occurs for 

„ „ f£ „ I = I* (22) 
c A 

for r = 10 9 cm and the wavelength A = 10 cm. 

The break down of the WKB method may be illustrated graphically using a different analogy. On a 
two-dimensional plane consider the surfaces of a constant phase of Hankel function Hy : 

<j> = arctan \ H , I = 0, 1, 2, ...v - 1 (23) 



Jy(Ar) v ) 

(Fig ||). This whole picture is rotating with a frequency lu. The points when the argument equals order, 
i.e., the "classical reflection point", may be regarded as light cylinder. For radii much larger than the "light 
cylinder" radius, r v = v '/A, the surface of a constant phase has a form of unwinding spiral with a wavelength 
\jv. Near the light cylinder the radial wave length becomes comparable to the light cylinder radius. 

From the quantum mechanical point of view, the kinetic energy, given by the derivative term in Eq 
. ( |TT| ) becomes zero at the light cylinder radius. This is the classical reflection point. For larger radii 
the kinetic energy is real and positive: the solution of the "Schrodinger equation" ( |i"i"| ) has a wavelikc 
structure. For smaller radii the kinetic energy is complex: the solution of the "Schrodinger equation" ( |TT| ) 
has exponentially decaying form, modified by the presence of the of the reflection point at r = 0. 



2.2. Wave-particle interaction 

A considerable simplification of the wave-particle interaction in cylindrical coordinates can be made 
when the interaction occurs in the WKB limit, or, equivalently, when we can use the expansion in tangent 
for the Bessel functions. From Eq. (|l^) we find a lower limit on the azimuthal wave number: 

v > (Ar) 2/3 w 10 5 (24) 

for the observed wavelength 10 cm and the curvature radius 10 9 cm. From Eqs. ( |22| ) and ( p4j ) it follows, 
that for a given radius and given wavelength cylindrical waves with azimuthal wave numbers in the range 
(Ar) 2 / 3 < (Ar) can be considered in the WKB limit, i.e. in the plane wave approximation. 
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It is possible to picture graphically the resonant and nonresonant interaction of a particle moving along 
the circular trajectory. The pattern in Fig. || may be thought of as rotating with a frequency Xc/v. A 
particle at a given distance r with a given velocity is rotating with an angular frequency f2 = v<p/r a . A 
resonant interaction occurs when these two frequencies are equal. In vacuum, 

r = « r v (\ - (25) 

where r v = v j\ is the light cylinder radius for the mode v. So, in vacuum, the resonance between a particle 
and a wave always happens inside the light cylinder. 



3. Synchrotron Emission Considered in Cylindrical Coordinates 

In this section we find electromagnetic fields from a particle executing helical motion in vacuum 
and rederive the conventional expressions for the cyclotron (or curvature) emissivities using cylindrical 
coordinates. The results of this section shows how the conventional expressions for the synchrotron 
emissivity in vacuum can be obtained using our approach. 

We first find the expansion of electromagnetic fields due to a particle executing helical motion in 
vacuum and in a homogeneous dielectric in terms of normal vector modes in cylindrical coordinates. This is 
done using dyadic Green's function that gives a response of a vector field to a vector source. We then show 
how conventional expressions for synchrotron or curvature emissivities can be obtained by two different 
methods: first by calculating an emissivity into a given cylindrical mode using an expansion of a current 
in terms of the normal vector modes and secondly by calculating Poyting flux through cylindrical surface. 
These calculations show how the resonant wave-particle interaction can be reformulated for cylindrical 
waves for a particular choice of helical motion . 



3.1. Dyadic Green's Function 

In order to find the electromagnetic fields produced by a given charge distribution in a medium is it 
necessary to know the Green's function for the vector d'Alambert's equation: 

□F(r,i)=Q(r,i) (26) 

The solution of the Eq. ( p6| ) may be written using dyadic Green's function Q(r, t, r', t') 

F{r, t)= J dv'dt'g{v, t, r', t') • Q(r', t') (27) 



In what follows we limit ourselves to the sources changing with time cx e tuJt . Then in vacuum 
d'Alambert's equation reduces to Helmholtz equation 

uj 2 

curl curl F(r, lo) -F(r, oj) = Q(r, ui) (28) 

c z 

where F is electric or magnetic field and Q is a source function. For a single particle source, which is a 
delta function in coordinates, Eq. j28|) becomes an equation for the Green's function. To find Green's 
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function we first have to find the solutions to the homogeneous equation with Q = 0. The solutions to the 
homogeneous vector wave equation can be obtained from the solution of the scalar wave equation 

uj 1 

curl curl f T ^ = (29) 

k x , A) = e^+^Z^Xr) (30) 
(Z v is any Bessel function) using the relations 

E (t) = curl(e x *), E {lt) = curl curl (e^) (31) 
Functions (|30|) arc cigcnfunction for the scalar wave equation. The eigenfunctions for the vector wave 



equation ( p9|) are (e.g. Eq. (6.5-6.6) of Tai 1994 ) 



L = igrad* = eW*.*) f Z ' v e r + -^e + +^Z„e x ^j (32) 
M = -j-curl^*) = e W^) (^Z u e r - Z'„e^ (33) 

N = 4curlcurl(e x vl/) = c l e K^ x ) (^ Z ' v e r - ^Z^ + Z v e x ) (34) 
uj\ uj \ A \ A r J 

We note that curlL = 0. L is an expansion in eigenfunctions of the longitudinal electric field from the point 
source and M and N are expansions in eigenfunctions of the transverse fields. 

The orthogonality properties of functions ( p2] - |34] ) are 

/ L*(r,A',fc>') \ / L(r,\,k x ,v) \ S(X-X') 
dr\ M*(r,X',k' x y) M(r,X,k x ,v) = (2tt) 2 (1 + S v , ) 1 ' S(k x - k' x )8 v - v . (35) 

J \ N*(r,A',£4X) / V N(r, X,k x ,v) J 

For the particle moving along trajectory r (t) the eigenfunction expansion of the scalar Green's 
function for the scalar wave equation 

curlcurl*-^* = -47r / dte iujt ^ - r °®h (x - x (t))8(<f> - (36) 
cr J r 

is 

ni \ u \ a xi o u \ J Ju{Xr)H [ u ] ' (Xr ) if r < r 
&(r,r ,X,v,k x ) =ind(uj- vll-k x v x ) < N „m,, \ . r ^ ( 37 ) 

[ J v (Xr Q )Hy(Xr ) if r > r 

This particular choice of Bessel and Hankel functions insures that the solution is regular at zero and 
corresponds to the outgoing wave at infinity. 



The corresponding dyadic Green's functions for the vector wave equation (|28j) can be found using the 
relations 

g B (r,r Q ) = curlXG(r,r ) (38) 
G E (r, r ) = (?+ ^ vv ) G ( r < r o) ( 39 ) 

where X is a unity matrix, G(r, ro) is the scalar Green's function and <?s(r,ro) and <?£j(r, ro) are magnetic 
and electric dyadic Green's functions. 
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The eigenfunction expansion of the electric dyadic Green's function is 
Qe{t, r , A, v, k x ) = in(2 - S Vt0 )5 (u> - k^v^ - k x u d ) 

L*(Ar)®L(Ar ) + M*(Ar)®M(Ar ) + N*(Ar)®N(Ar ) I (40) 



which, using relations ( p9| ) gives 

Qs(r, r , A, v, k x ) = in(2 - 8 vfi ) 

(#«'(Ar)e r - ^F«(Ar) e<# - ^fl^Ar)^ ® K(Ar )e r + ^-J„(Ar o )e + ^ J v (Xr )e x j 
+ C -£- (^y^'^'K - ^^ 1} (Ar)e, + ff«(Ar)e x ) 
(g) f i-^-J' u (\r )e r - j^- J^(Ar o )e0 + J i ,(Ar )e ;E J 

+ (^ff«(Ar)e r - £T«'(Ar)e ) ® J,(Ar )e r - 4(Ar )e^ J (41) 

The first term here is just the eigenfunction expansion of the Green's dyadic for the point source. 

The magnetic Green's function (the Fourier transform in 4>, x and t) for the vector wave equation is 

r , A, v, k x ) = in(2 - 5 U>0 ) | N*(Ar) ® M(Ar ) + M*(Ar) ® N(Aro) | (42) 

which, using relations (p9[) gives 



cA 

GB{r,r Q ,\,v,k x ) = iir(2 - <5„ ) — 



(-ij-HW(\r)er - H^'(Xr)e^j <g> ^^(Ar )e r - J,(Ar o )e + J,(Ar )e x 
+ f-i^ffM'^K - ^fl-W(Ar)e^ + ffW(Ar)0 ® ^ J,(Ar )e r - ^(Ar )e,, ) ) (4;-!) 



The solution to the inhomogeneous vector wave equation ( p8[ ) is given by the integral with a kernel 
given by the dyadic Green's function 



F(r, A, v, k x ) — J r'dr'g(r,r',\,v,k x yQ(r',\,v,k x ,co) (44) 
where Q(r', A, v, k x ) is a Fourier transform of Q(r). 
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3.2. Synchrotron Emissivity into Cylindrical Mode 

In this subsection we show how conventional emissivities into plane waves can be modified to obtain 
similar relations for the emissivities into cylindrical modes. To simplify the analysis we consider in detail 
only emissivity of a charge moving along a spiral trajectory. An important factor that simplifies our 
calculations and allows one to obtain the emissivity in a simple closed form is that for helical motion a 
particle always stays at the same radial coordinates. If the particles were not on the ground gyrational 
level, then in calculating the emissivity it would be necessary to average the gyrational motion of particles 
taking into account the radial dependence of the vector normal modes. 

Following the conventional approach in the theory of electromagnetic wave-particle interaction (e.g, 



Melrose 197q) , we can identify the source of energy in the wave with the the work done by the extraneous 
current: 

P = J drdtE(r,t) ■ j(r,i) (45) 
Using the power theorem for the Fourier transform, i.e., 

/dk 
^f(k)g(k) (46) 

(where f(k) and g(k) are Fourier transforms of f(t) and g(t)) we can rewrite Eq. ( ^H|) for the Fourier 
transforms of E and j in time, <j> and x Q 



w) (47) 



Electric field in Eq. (E^) may be expressed in terms of the dyadic Green's function and the current (see Eq. 
([44|)). For a particular case of a particle executing helical motion with velocity v = {0,v<p,v x }, the Fourier 
transform of the current density j(r) = gv<5(r — r ) is 

j(r, A, v, k x , u) = / dxe~ lk * x [ dte wt [ # e -^j(r) = qv5(cj -vQ- k x v x ) 5{r ' T °\ (48) 



so that the electric field is given by (v ^ 0) 

E(r, A, v, k x ,u) = ^ J2 Ni*(Ar)(Ni(Ar ) • v)5(uj - vtl - k x v x ) (49) 

i 

where Ni = L, N, M Inserting Eq. ps| ) in Eq. ( [47| ) we find the radiated energy 

p= 2^^ |Ni(Aro) v|2 (5Q) 



which gives the cyclotron emissivity 



(51) 



which is exactly the cyclotron emissivity per unit interval d<p dtu dk x . 



1 it is also possible to make also expansion of the vector fields E and j in terms of the normal modes L, N, M and use later 
the dyadic Green's function G(X,u,k x ,u>) 
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3.3. Alternative Derivation of Cyclotron Emissivity 

We can use an alternative method to derive the cyclotron emissivity into cylindrical modes by finding 
the electromagnetic fields produced by the given charge distribution and integrate the Pointing vector over 
the cylindrical surface. Using Eqns ([ll]), ( |44| ) and jig) we find the transverse component of the electric 
field: 

E x = nq(2 - <y„, )~J(w -vSl- k x v x ) (-J' v (\r )v^ (-i^if«(Ar)e r - '(Ar)e ) 

(-i^Wer - gijW(Ar)e, +i/( 1 )(Ar)e^) (52) 

In the wave zone, Ar ^> 1, and for v 3> 1 this can be simplified 



, 2 a; , A „ / T/ . . . c 2 A / fc^j/O \ , [ k x 

Ej_ = «97ta/ ^:~ e I J v (Xr )v^ - — — I ^ + ^ J ^( Ar o) I *~ e '' + e * 

<5(ci> — t^ri — k x v x ) (53) 

Similarly, we find for the magnetic field 

B = inq(2 - £„ >0 )A (^(-z^tf «(Ar)e r - ff« '(Ar)e ) (~^% + «*J ^(Ar ) 
- (-^^'(Ar^ - |^ff«(Ar)e + ff«(Ar)e^ J' v {\r Q )v^ 8{lo - vil - k x v x ) (54) 
which in the limit Ar 3> 1 and v 3> 1 gives 

B = iqn\j -^-Ae lAr <5(u -vVt - k x v x ) (iJ v (\r ) ( h^l _ v \ e< ^ + (~^-e r + e x ^j J' v (Xr )v^\ (55) 

The Pointing flux, S = cE* x B/ (An) for r — * oo is then 

S = ^-^ {^J'vv\ + (^-^- - -j^j ^(Ar )^ ^e r + ^-e^ j <5(w - ^ - k x v x )dldwdk x (56) 

where = rd<^> is a unit arc length of a cylinder. 

We are interested only in the radial component of the flux. Equating the Pointing flux to the emissivity 



we find the cyclotron emissivity (51) 



4. Waves in an Isotropic Dielectric 

We expect that the dielectric properties of a medium, i.e., collective effects, will play an important in 
defining the properties of the wave-particle interaction. So, as a next step, we consider in brief the waves in 
an isotropic, frequency dispersive medium. In such a medium the wave equation is 

<9 2 D 

curl curl E - — — = (57) 
at 2 - 
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In an isotropic, dispersive medium the relation between the electric induction D and and electric field may 
be written as 

D(w) = e(w)E(w) (58) 
From Eqs. (g) and (^8|) it follows, that in a dielectric, the solutions of the wave equation are 

Ef)(r,fe,^) = ( l ^j'Xxr)e r - ^J u (\ r )e* + J„(Ar)e x ) (59) 
V \ 2 X 2 r / 

BC'M.V) = f ^e r + ! j:(Ar)e,) (60) 

A 2 c V r / 



E«(r 1 fe ) i/,u)=S( < ) ( -J„(Ar)e r - J^) (61) 



BM(r, a;) =.EW ( ^(Ar)e r + ^J„(Ar)e - J„(Ar)e x ] (62) 



where A = \J uj 2 /c 2 e — k 2 . 



4.1. Cherenkov Emission in Cylindrical coordinates 

As another didactic example we illustrate the features of the fields in the case of Cherenkov type 
emission taking as an example scalar waves. The solution for the vector waves can then be obtained using 
the general relations between the solutions of the scalar and vector wave equations. 

In cylindrical coordinates the equation for the eigenfunction expansion of the scalar Green's function is 
given by Eq. ( |37| ) with vq = 0. The only remaining term is 

G((r, r , A, v = 0, k x ) = iir5(u> - k x v x )H (Xr) (63) 

Using the argument of the (5-function we find 



\=\j—-k 2 x = k x ^i~l (64) 

where (3 X — v x /c. For e(3 2 < 1 the argument of the Hankel function is complex, so that the fields decay 
exponentially at large distances. For eft 2 > 1 the argument of the Hankel function is real, which corresponds 
to outgoing waves. For e/3 2 = 1 the Green's function ( j63| ) has a discontinuity, this is the shock front 
corresponding to the Cherenkov cone. 



5. Dispersion Relation in an Infinitely Strong Curved Magnetic Field. Vlasov Approach. 

The next problem that we will consider is a dielectric tensor in cylindrical coordinates for plasma in 
the infinitely strong circular magnetic field. In this particular case, it is possible to relate the electric field 
and the electric induction in cylindrical coordinates through a conventional dielectric tensor and not the 
dielectric tensor operator, in other words the dielectric tensor becomes local. 
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5.1. One Dimensional Plasma in Inhomogeneous Magnetic Field 

The equations governing the electrodynamics of a plasma in a strong , weakly inhomogeneous magnetic 
field are the wave equation 

> <9 2 E(r,t) 47r<9j(r,f) 
curl curl E(r, t) + ^2 = ~ dt (65) 

and collisionlcss Vlasov equation 

Bf(r, t, P) +(v .|. )/(r , t)P) + ,f E+ i vx(B + B )^ = (66) 



dt dr ' ' \ c J dp 

where B° is the external magnetic field, /(r, t, p) is a particle distribution function, j(r, t) is a current, q is 
a charge of particles. 

We represent the distribution function as a sum of the unperturbed distribution and fluctuating part: 
/(r,i) = /(°)(r) + /W(r,t). We assume that the unperturbed distribution function is stationary, i.e., 
df^ /dt = 0. A real plasma present on the open field lines of the pulsar magnetosphere may not satisfy this 
condition. It can be hydrodynamically unstable with a growth rate much smaller than the dynamical time 
(which is equal to the pulsar period). But since in our model problem, the plasma circulates infinitely long 
along the circular magnetic field, this condition is necessary to impose. Then the unperturbed distribution 
function satisfied the equation 

In other words, /^°^ must be a function of the integrals of motion. The integrals of motion for a particle 



:/ m + G-*°)tr = ° m 



in azimuthally symmetrical static magnetic field have been considered by Mikhailovski 1992. The relevant 
constants of motion are the radial coordinate of the gyration center 

r g = r + ^, (68) 
parallel velocity of the gyration center V^, which, up to terms linear in the small parameter of the adiabatic 



approximation coincides with the parallel velocity of the particle V$ = and the quantity Vj_ ( Mikhailovski 



1992, Eq. (15.11)), which, in turn, up to terms linear in the same small parameter coincides with the 
gyrational part of the particle velocity V± = v± (see below) . The most general form of the distribution 
function is then 

fW = F(Vx,Vi,r g ) (69) 

where F is an arbitrary function. We also note, that since the distribution function is independent of x and 
4> by our choice and cannot depend on the gyrational phase of the particle, the first term in Eq. ( |67| ) is zero. 



Following Mikhailovski 1992, we find, that the inhomogeneity of the external magnetic field results in 



an equilibrium electric current (diamagnetic current) given by ( [Mikhailovski 1992| ) 



B V dr 

where P± and P» are transverse and parallel plasma pressures, defined by 



p \\ = E / W/ ( ° )dv (71) 
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Here v\ = v 2 + (v x — u<i) 2 , Ud is a drift velocity. 

Using Maxwell equations and Eq. ( f70| ) the condition of the transverse equilibrium becomes 
( [Mikhailovski 1992| ) 

d ( B 2 \ B 2 P±-P\ 



Or \ 07T / 4-7T T 

In the case P± — the diamagnetic current ( |70|) may be written as 

j x = gnud, u d = (73) 

Lu B r 

where Ud is the drift velocity directed along the binormal to the magnetic field. We note here, that there is 
an alternative way of deriving Eq. ([73|) from the single particle equation of motion using averaging over fast 
rotations around magnetic field (e.g ). After averaging over fast rotation the resulting equation of motion 
for the location of the gyration center contains inertial forces due to the inhomogeneity of the magnetic 
field. These inertial forces result in a drift motion which produces diamagnetic current. 

Returning to the Vlasov equation, we make a Fourier transform in time, 4> and x (oc exp{i(y(j)-\-k x x— u>t)) 
and find an equation for /W (^ m, k x ,u>): 

i(-ui + + k x v x + v r — )/ (1) (r, m, k x ,oj,p) + q ( E(r,m, k x ,w) + -v x B(r,m,k x ,w) ) — — = (74) 



or \ c J op 



5.2. Infinitely Strong Magnetic Field 

Equation (|74j) is a differential equation for f^ (r, to, k x , ui, p). This reflects the fact, that in cylindrical 
coordinates (any curvilinear coordinates) the relation between electric field and electric induction involves 
dielectric tensor operate^ The case of infinitely strong magnetic field allows a considerable simplification: 
in this case the operator relation is reduced to algebraic, i.e. , to the conventional dielectric tensor. 

In the case of infinitely strong magnetic field the velocities across the field are zero: v x ,v r = and Eq. 
( f74| ) can be solved for (r, to, k x ,oj): 



/(l>(7,TO,^,p) = ^^^ 

C(UJ — His) 



(75) 



The corresponding current density is 



3<f> 

where 



fV{r,m,k M4> )vtdpt = ^ J dp. ^-^^L- E, (76) 



is a one dimensional distribution function. 



2 Mathematically, the difference is that dielectric tensor acts in a space tangent to the vector field at some point, while 
dielectric tensor operator acts on a vector field itself. 
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We can now introduce a dielectric tensor 

/ 1 

e l3 (r,v,k x ,Lu) 1-K | (7S) 

V o o i 

where 

K = / p-^ = / dp^ ^ — 75 79 

Note, that both ji(r, v, k x , u>) and Ej(r,is,k x ,Lu) are conventional Fourier transforms, so that the hermitian, 
e^(r,v,k x ,Lu), and antihermitians, ef-(r,v,k x ,Lu), parts of the dielectric tensor ( f78| ) satisfies the usual 
conditions 

eij(r,v,kx,w) = e"(r,-v,-k x ,-oj) 
£?j(r,v,k x ,uj) = -ef!j(r,-i>,-k x ,-w) 

eij(r,-v,-k x ,-w) =ey(r,-i/,-fc !B ,-w) (80) 

The wave equation now reads 

'" 3 (r^) + i*^, + f ^ - f 4 - ^ = (81) 



dr 

d ( 1 d , „ . \ . d ( E r \ k x v „ / u; 



,2 



19 ^V-l-w-f^--)^ " 



r dr \ dr J r dr 

To illustrate the solutions of Eqs. (^l| - |S3|) we show how, in the WKB approximation, the solutions 
can be obtained for a given dependence K(r) in the case k x = 0. Then we find exact solutions for k x = 
and constant (independent of r) K. Unfortunately, there seems to be no simple way of solving Eqs. (|3|) in 
known functions for k x ^ 0. 



5.3. Propagation with k x — 

For k x = the equation for E x separates from the equations for E r and E^. From Eq. ( |S3| ) we find 

E (lt) (r,4) =exp{i(v(f>-ujt)}J v (^) e x (84) 

B (lt) (r, t) = - exp{z(z/0 - ut)}^ (^J v (^) e r + i J' u (^) e ) (85) 



Using Eqs ( |84| ) and (85) it is possible to show that the following relation holds: 

[l-K)vE <j> -i^-{rEr) = Q 
or 

which can be derived directly from the relation divD = 0, where Di — tijEj is the electric induction. 
Using Eq. <M) we find 



(86) 



—E - 

dr 2 r r 



r dK\ dE r 



1 — K dr I dr 



1 1 
(1-K ) + - + —— 

r z r[l — 



dK 



(1 - K) dr 



= 



(87) 
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Equation (pq) can be solved using WKB method for a given dependence K(r). If K{r) is a constant 
(as a functions of r), then it is possible to find exact solutions: 

E< 2 )(r, t) = exp{i(v<f> - cut)} (^J v ^k (v^ 1 ^) &r ~ ^VT^ ^^~^) e *) 
BW(r, t) = exp^W - urf)}^J (^(1 - K) + J v ^ (^VT~??) (89) 



5.4. Curvature Emission in Vacuum and in the Isotropic Medium 



In this section we'll show how the conventional emissivities for curvature emission can be obtained 
in this approach. We will calculate the single particle emissivity as the energy dissipated by the optically 
active medium. This can be done by relating in Eq. (ET]) the local current j(r, v 1 k xi ui) and the electric field 
through the antihermitian part of the dielectric tensor ([78j) : 



ji{r,v,k x ,uj) = -i—e^Ejir, v,k x ,u) 

4-7T J r 

To procede further we define the wave amplitute, E-^. , as 

E(r, v, k x , u) = 2tt J2 ^N i N i (r, V) k x ,X)S{X - y/rf/c 2 - kj) 



(90) 



(91) 



Then the energy density in the wave W is (V — tt1 2 H) 



i r T/2 i r H/2 . ^ i r l , iE(r,oi 2 



W — lim lim lim — . 

T — >ool — >oo H — >oo 1 ./_y/2 J — Hjl Jo 



J dk x duj 



f rdr 1 



lim , n . „ / dk x duj 
Moo^ (2?r) 3 ^ J 



(II— I il.r 

E(r, v, k x ,uj)\ 2 



fy— / rdr- 



8vr 



8tt 



— |Ni| 2 [5(A- Va;Vc2-fc2) 

7Tt z 



(92) 



Using the relations 



lim / rdrJ 2 {Xr) — lim — 

l — >oo ^ / — >oo 2 



J^(Ar) + 1 



\ 2 f 



J 2 (Ar) 



7rA 



lim y[5(A - V^T^fcf)] 2 = ^(A - ^ 2 /c 2 - M) 

l— >oo I Z7T 



(93) 



We find 



^ = E / dk x dJ-^6(\ - V^Jci—k 2 ) = 7A3 E / ^Wk, 



(94) 



where we identified an energy in a vector mode Nj as 



(95) 
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Dcfining the energy gain of the wave as —TiWi we find using Eq. (47) 

-Y^J^yll f dk * dLuT * W * = -*^J^3 E J dk x duj\E N; | 2 J rdrN-*(Xr) ■ e A ■ N ; (Ar) ^ ~ ? ' o) (96) 
so that 

r i --N i *(Ar )-e A -N i (Ar ) (97) 

7T 

This expression provides a growth rate for a vector cylindrical mode Ni = N, M due to a current circulating 
at a radius tq. 

To find the single particle emissivity we assume that the growth of the wave is due to the resonant 
interaction with plasma particles which emit and absorb wave quanta. We then can identify the single 
particle probability of emission w{p,p, k x , A) as 

/Q P (p) 
dp c / ) w(p t p,v,k x ,X)hk^—^- — (98) 

Using the antihermitian part of the dielectric tensor 

e <p<t, = -1 / dp^v^— S (u - vVt) (99) 



iira 2 vl / u 2 k 2 r 2 

w(pt, u, k x , A) - — J' 2 (\ r ) + ^-Jl{\r) ) 8 (u - rtl) (100) 



We find 

47T0 2 ^ / ,„ 
.nCr., u \\ = - t I _T 2 l\ r \ 4- 

\ 2 1 2 

and emissivity r](v, k x , uj) 



This is emissivity per unit range dk x dujd(f>. 



M == E ( J '» 2 ( Xr °) + x^ J " {Xr) ) 6{UJ - M) (101) 



Effectively, we calculated the induced energy loss of a particle due to the interaction with a wave and 
then, using Einstein relations between the induced and spontaneous emission, we find the spontaneous 
emissivity. 

Introducing notations k x — a>/csin#, dk x — lo jc cos 8d9, A = w/ccosO and = Qr/c the expression 
for the curvature emissivity may be reduced to [] 

V (u>) = [J' v 2 (p^cos0) + ^g^J 2 (^cosc?) J Siu-vp+c/r) (102) 

This is emissivity per unit solid angle dO. This is exactly, the expression for the single particle curvature 
emissivity in vacuum ( Melrose 1986 , Eq. (13.62)). We stress, that Eq. (102) was obtained in cylindrical 
coordinates. It gives a single particle emissivity per unit frequency, which is, of course, independent of the 
coordinate system used. 

As the waves propagate outwards from the emitting region, the approximation of an infinitely long 
cylinder will cease to be true. Then, at some transition region the cylindrical waves will become spherical. 



3 Since angle 9 is measured from the plane k x = the unit solid angle is dVl = d(f> cos ddd 
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It should be possible to consider the details of such transition using Kirchoff integrals over a remote surface 
where dispersive properties of a medium are not important. Here we only note that the parameter 6 
introduced above corresponds in a vacuum case to the polar angle in spherical coordinates. This provides a 
simple formal relation between the spherical and cylindrical coordinates. 

In our approach, which treats the wave-particle interaction as resonant, the impossibility of the 



amplified curvature emission in vacuum, first derived by Blandford 1975, follows from the fact that in 
vacuum there are no subluminous waves, so that no waves have a phase velocity that could fall into the 
region where the derivative df(p < p)/dp^ > is positive. 

Similarly to vacuum case, we can obtain a curvature emissivity in an isotropic medium: 

v(to) = ^ ^4 2 (^/3 cos0) + ^/j2(^/3 cos0)^ 8(w-v0tc/r) (103) 

(emissivity per unit solid angle dQ). 



6. Airy function approximation 



In this section we consider the Airy function approximation to the curvature emissivities of 
ultrarelativistic particles in a medium, which is a common way of approximating the cyclotron and 
curvature emissivity of ultrarelativistic particles in vacuum (e.g., Melrose 1978|). We show, that for a 



particles moving with the velocity larger than the speed of light in a medium which has e > 1 We assume 
that a particle is moving through a dielectric with the permittivity e > 1 a qualitatively different expansion 
of the Airy function should be used. 

In the transition region, when the argument of the Bessel functions is large and close to their orders it 
is possible to use the Airy function approximation 



1/3 



J v (v + zv 1 '*) « \ ±\ Ail 2 1 :! : 
where Ai (x) is Airy function 



j2/3 
3;/ 



2 2 / 3 y/I 



VW/a Kl / 3 



3/2 



if z < 



(104) 



In vacuum, for a resonant wave-particle interaction the argument of the Bessel functions is always 
smaller than the order: Xr — v(3,p cosd < v. In a medium, the argument of the Bessel functions can become 
large than the order. In a dielectric, the argument of the Bessel functions is Xr = yJoj 2 e/c 2 — k 2 r. If we 
introduce a notation k x — uy^e/csuiO (this is a definition of angle 9, the condition of regularity at infinity, 
A 2 > insures that sin9 < 1), then we have Xr = v^fift^ cos 9, which can be larger than v if > l/y/e - 
for the superluminal motion of a particle. Using these notations we find 



z = (v^oos* - « (6- -L - ^) v 2 ^ 



(105) 



for e = 1 + 25, 6 <C 1, 7 S> 1 and <C 1. It is clear from Eq. (105), that in vacuum z < and in a medium 
z becomes positive for superluminal particles. 



Following the discussion of Section 2.1.2, we can identify the light cylinder radius for the mode v as 



r„ = v / X. We can argue, that in the case z > the resonant interaction of a particle with a wave occurs 
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outside the light cylinder 

r ° = ^{ 1 + £s) ( 106 ) 

A transition through point z — ("light cylinder" or "classical reflection point" ) is nontrivial. It 



resembles phase transition (Bchwinger et al. 1976) in a sense that correlation length for thermal or quantum 
fluctuations is very large near the transition. The physical conditions beyond and above the transition point 
are essentially different. 

For superluminal motion the collective effects of the medium play important role. In this case, the 
corresponding emission process can be called Cherenkov-curvature emission, stressing the fact that both 
inhomogeneous magnetic field and a medium are important for the emission. In the vacuum limit, n — ► 1, 
Cherenkov-curvature emission emission reduces to the conventional curvature emission. Conventional 
Cherenkov radiation may be obtained in the limit r — > oo after integration over v (sec [Bchwingcr ct al. 1976 
for the corresponding transition for the cyclotron-Cherenkov radiation). 



The emissivity for the Cherenkov-curvature process is 

*0 - ^Sf 1 (^W) + (!)""•(-»■/••>)<<.,-.*> (107, 
For z < this reduces to the conventional representation for synchrotron emission in terms of 



MacDonald functions Ki/ 3 . For z > Eq. (107) gives 



2 1 l z JlJ 1 q 2 v 2 , /tan 2 6 9 2 2 / 3 z 2 a\ 

= w/3 OM0 + J-irfQ) + — (-M0 - J-ys(0) j a (« - **) (ios) 



where we used 

2/3 



J' u (v + zu 1 ' 3 ) « (?) At' (-2^z) 
Ai'{-z) = - Z - (j_ 2/3 (-2V3*) - J 2/3 (-2 1 / 3 z)) (109) 



Emissivity (108) may be simplified in the case ^ 3> 1. Then we can use the asymptotic expansion for 
Airy functions 

Ai(-z) 



,1/4 
J* 

to find 



(r 3/2 


+ 


7T 






4 


jz 3 / 2 + 


7T 


) 




4 





Ai ' ( _ z ) ra £_ COB (^3/2 + « ) (uo) 



V2V~eu 2 q 2 PlVi ( 3ft , ... , sin 2 (£ + V4)*/ 2/3 tan 2 fl\ , n . 



where £ = (2z) 3 / 2 /3 



For (5 ^> I/7 2 , and 5 ^> 9 2 the condition z » 1 implies Sv 2 / 3 ^> 1. Then Eq. (Ill) can be further 
simplified: 

V2^V5co 2 q 2 p 2 . (2 z l 2 v tt\ 
r?M « V ^ cos 2 (^^-^ 3/2 + 4 j ^ - "») (112) 
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In this context we note, that the total spectral power for the curvature emission in a medium, viz 



V (oj) = ^ 2n 2 (3l J' 2v {2vnM - (1 - ti 2 /3 2 ) 



(113) 



for the case of superluminal motion, n(3^ > 1, can be reduced to the explicitly Cherenkov-type emission 
form: 

1 



with A(z) of the order of unity for z > 1 ( Schwinger et al. 1976| ) 



A(z) 



(114) 



7. Resonant Electromagnetic Waves in the Asymptotic Regime z ^S> 1 



We wish to simplify the relations (|59|-|62[) in the limit z 1. Instead of standing waves, described by 
Bessel J functions, we consider propagating waves, described by Hankel and functions. For the 
outgoing waves, corresponding to ifW , we find in the limit z 1: 



H {1)> 



2 l/4 



^l/3 z l/4 



_. 2 3/V/ 4 _^,3 /2+lf 



Next we define the radial wave number 



_ -e 3^ 

7TI/ 2 / 3 

.dlnifW 
<9r 



Using Eqs (115) and (116) we find 



,1/3 



25A 



(115) 
(116) 
(117) 



Introducing a notation — v jr the waves in the isotropic dielectric in the limit z 1 may be written 

(118) 
(119) 

E w (k,w) =E {t) (fc e r -fc r e ) 



E( !i )(k lW )=# f_^e r -^e + , 
V A 2 A 2 



B<">(k,w) = (fc e r - k re4> ) 

A 2 c 



B« (k, oj) = £?W [ ^e r + 



cA 2 



(120) 
(121) 



with k = {fc r , fc^, k x } and a phase dependence exp{— i(u>t — k r r — k^cj) — k x x}. 
With these notations the dispersion relation is 



+ k r + fcj. — 



2 _ W 6 



(122) 



Relations (118- 121) look similar to the relations (|16[) , though they were obtained in quite a different 
manner. Expansions (|16D are valid for Xr — v > v with the radial wave vector defined by ( |l4|). E xpansions 
(118- 121) are valid for Arwy and £ ^> 1 (the radial wave vector in this case is defined by ( |117| )). 
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The importance of these results is that in these two limits, the electromagnetic wave in cylindrical 
coordinates look like plane waves. This is a considerable simplification. It allows one to implement a well 
developed technique of Fourier transforms in considering the wave propagation. This representation of 
the electromagnetic fields requires (i) the presence of a medium with the index of refraction n > 1, (ii) 
superluminal motion of the resonant particle S > l/(27 2 ), (iii) large harmonic number v > <5~ 3 / 2 . 

Another important feature is that in this limit there is an additional freedom in the choice of wave 
polarizations in an isotropic dielectric. For example, we can introduce the polarization vector corresponding 



to the waves (|118| - 121 ): 



eW = i 
A 



— k r , 0} = {sin-0, — cos-0, 0} 
{— cos ip sin 6, — sin tp sin 9, cos 9} 



which satisfy relations 



e^ ^ — e^ x b x 
e<'*> = e k x e W 



(123) 



(124) 



The special role played by the in defining the polarizations of the normal modes comes from the 
fact that in the isotropic homogeneous medium the solutions of the vector wave equations may be chosen 
to be tangential to the coordinate surfaces e x . In the nonisotropic medium this is not true (see Sec. 

Alternatively, in the limit z> 1 we can chose the following polarizations 

1 



,(«) 



k± 
1 

kki 



{ k x , 0, kr\ 
'{k(pk ri kj_ : k x k(p} 



(125) 



(see Fig. ||). This particular choice of polarizations has an advantage that it may be related to the 
polarizations in the straight field line geometry. 



8. Response Tensor for a One Dimensional Plasma in a Curved Magnetic Field 

We will calculate the response tensor of a one dimensional plasma in cylindrical coordinates using the 
analogy of the forward-scattering method (see e.g. Melrose 1986 ) adopted to cylindrical coordinates. We 
represent the current as a sum of currents due to the each single particle moving along its trajectory: 



j(r,f) = J dr°dpi sp (r,t)f(p : r°(t)) 

We expand the single particle current (denoted by the subscript sp) 

j sp (v 1 t)=qrS(r-r°(t)) 
in Fourier amplitudes in time, k x and <j> and in Hankel amplitudes in r 

j sp (r,m,k x ,uj) — q J dtexp{-iu;t} J dxexp{ik x x} J d(j}exp{iv(j)} j sp (r,t) 
dtexp{i(v<P°(t) + k x x a {t) - ojt)}r / ^^(^^(^{t)) 



(126) 



(127) 



(128) 
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where we used a Hankcl transform of the delta function Q 



5{r-r°) 



We expand the orbit of the particle in powers of wave amplitudes: 

r(t) =r (i)+rW(i) 

The first order Fourier transform of the single particle current is then 

jW(r,m,^,w) = q J dtexp{i(v<f>°(t) + k x x {t) - wt} J ^d£_J u (£_r) 

The orbit of a particle is found by solving the equation of motion 

dp 



(129) 



(130) 



dt 



F°(t,r,v) + FW(i,r,v) 



(131) 



(132) 



where F°(t, r, v) is a force acting on a particle when no wave is present and pW(t, r, v) i.e. a force acting 
on a particle due to the presence of a waves. Expanding Eq. (|132|) in powers of waves amplitudes we find 



dp°(t) 

dt 
rfpW(f) 



= F°(t,r°,v ) 



(133) 
(134) 



The force F^ (t, r°, v°) can be expanded in Fourier amplitudes 

FW(t,r°,v°)=5Z / ^'*4exp{-i(w , t-fc> (t)-i/>°(t))F( 1 )(a;,m ) A; a: ,r (f),t) (135) 



Equation (134) can be solved for the first order velocity perturbation rW(t) and first order trajectory 
perturbations x^(t), (1) (f) and r^(t): 



( WW \ 

0(D (t) 
V r(D(i) / 



J2 / dtyd*4 exp{-i(o/t - k' x x°(t) - v'(jP(t)) 



I V(uj',i/,k' x ,r°(t),t) \ 
x^{u',v',k' x ,r Q {t),t) 



(136) 



where tilde denotes Fourier transforms. 

The first order single particle current then becomes 

jW( r , m ,t,,) = g ^ / duj'dk' x I dfe^tt"-" ■)t-{k x -K) x \t)-(u-u ')^(t)))} I Wy ^ r) 



VJ, 



•„(£r°(t)) +r°(t) + *fc,£«) J„(£r°(t)) + g^M f (D 



(137) 



A Hankel transform is denned as -F(x) = J„(£o:)£g(£ ydyJv(£y)F(y) 
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Using this single particle current we can calculate the current ( J126| ). In the general case, current (137) 
may be related to the electric field of the perturbing wave through a generalized dielectric tensor: 

j(r, m, k x ,uj) = £ (r, m, k x , w) ■ E(r, m, fc x , w) (138) 

where £ is a dielectric tensor operator which involves a Hankel transform of the external current and partial 
derivatives with respect to r. This is different from Cartesian coordinates, where the electric induction is 
related to the electric field through a dielectric tensor 5 . To simplify the calculations we assume that the 
radial dependence of both the Hankel transform of the external current and of the perturbing electric field 
can be approximated by the plane wave form. For the nonresonant term in the dielectric tensor operator 
this is justified in the WKB limit, while for the resonant terms this is justified for subluminous waves in the 
limit 8 2/ 2 / 3 > 1. 

Thus, in the Hankel transform of the external current and in the derivatives of the perturbing electric 
field, we can identify 

d d 

g^^ik ro , fr^ikr (139) 

In this approximation the Hankel transform of the external current will reduce to the Fourier transform in 
r and the derivatives of the perturbing electric field are replaced by ik r . The dielectric tensor operator £ 
then becomes a conventional tensor. 

The integration of the induced current over d<fi°dx° then gives delta-functions 8(v — z/) 5(k x — k' x ) that 
are subsequently removed by the corresponding integrations. For tq = const, the time and du' integrations 
insure that only secular terms are retained. The integration over d£ and dr° are removed using relation 



8.1. Perturbed Trajectory 

The equations of motions for a particle in a circular magnetic field when the Larmor radius is much 
smaller than the radius of curvature are the following: 

^ = e£Ur°) + q - ((i? (r°) + B )v x - ^(r )) (140) 

^ = eE x (r°) + ± (B r (i>0 " v r (B^(r°) + B )) (141) 

^ = e£ (r°) + q - (B x (r°)v r - v x B r (r )) (142) 

We solve Eqs ( |14C| - |142| ) by expanding in powers of the wave amplitudes. We assume 
that initially the particles are in the ground gyration state. Then in the zeroth order we find: 

v 2 

v® = 0, <f> = Clt = v^t/r, v x — Ud, with Ud — Rb * b ■ The first order in wave amplitudes gives 

(£ r(r 0) _ V -l Bx (r°) + ^B,(r )) + (143) 

V C C / T 



dvf- q 
dt "/rrieC 



5 Mathematically, the difference is that dielectric tensor acts in a space tangent to the vector field at some point, while 
dielectric tensor operator acts on a tensor filed itself. 
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dv. 



(i) 



dt "/m e c 



c c 2 J 7 



dv 



(i) 



dt "/m e c 



' { 1 E^r ) - —B r (r°) - M v 4> E 4>( r °) + u d E x (r )) 



Expanding Eqs (142 - 145) in Fourier amplitudes we find the solutions: 

B^Ud — BxVij) 



TO e c7(cj^/7 2 — f2 02 ) 



in ( E, 



UJ B I , V^Br UdivjEj, + U d E x ) 



ul 1 ) = 



m e cy(uj B /"f 2 — fl°< 2 ) 

I EL 



C 

m° | e x 



V^Br Uddv^E^ + u d E x ) 



7 



B^Ud — B x v§ 



* m e fi°7 
where f2° = uj — k^v^ — k x Ud, k^ — m/r. 

The first order variations in trajectory are 

,(!) 



, r u d D \ v^v^E^ + UdE x ) 

— \&<b £>r\ H o 



,(!) 



In cylindrical coordinates, the relation between the electric and magnetic field is 



B r = 



•E, 



k<hC 



E , 



UJ UJ 

k x c .cdE x 

Ba = E r — i — — 

uj war 

B x = J^E r +i — —(rE< b ) 
ui uir or 



(144) 
(145) 



(146) 

(147) 
(148) 



(149) 



(150) 



8.2. Simplified Response Tensor 



Using Eqns (146 - 15C) we can find the dielectric tensor operator £. In Sections 2.1.1 and 2.1.2 we 
discussed when the operator relations between electric field and electric displacement can be simplified to 
algebraic relations for the nonresonant waves, and in Section [7] we found a limiting case z» 1 when this 
can be done for the resonant waves. In these limits it is possible to change all the radial derivative d r — ► ik r . 
The relations between magnetic and electric field then simplify 



B,< 



B, 



k x c k r c 

&r H &x 

UJ UJ 

k^c k r c 

&r &tb 

UJ UJ 



(151) 



The resulting dielectric tensor is still complicated. It can be simplified according to the following 
procedure: (i) for the nonresonant parts the drift velocity is small and can be neglected, (ii) for Cherenkov- 
type resonances (resonances that do not involve u>b) we retain all the terms, (iii) for cyclotron- type 
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resonances (resonances that do involve u>b) we assume that k x /k,p w u^/c, Ud/c^> I/7 and expand in 
small drift velocity retaining only the first terms in Ud and k x . Implicit in this expansion in small Ud is the 
assumption that 7U 2 3> u> 2 B u 2 d /c 2 . In this limit, in the cyclotron-type resonant terms involving drift velocity 
we can approximate v^, ss c. The dielectric tensor is then 

dp<t> fa k l U l ~\ 



= 1 ~ \ E ~S j ~ ~ kxUa ^ A «f a 



dp<f, fa | 



_ w p« /* dp* (k 2 x + k 2 )v 2 1 w 2 Q /■ ^ 

^ y 7 /a 24- ^ y 7 fcrUaA « /a 

= "5 E ^ / ^ ((" - M*)^ + <*r«*£) /« 
a J 7 

= J E ~¥ / ^ (( W - fc <A^) A a - **r«aA+) / Q 
a ^ 

\- Wp Q /" dp* (/c 2 + A: 2 ) 

-E w p« / 



/a k^u a j k x 1l a ^ 



7 ^S 2 w 
\Y^^j d -fv,{k x At- lkr A-)f a 

dp4> fa k^u a ( k x u a \ f luv^ 



-E"& J 



7 f2° 2 w \ u J \ k<j,c 



a J ' 

= - i E I (^1 + *k x A-) fa (152) 

Here 

= ( — + tt^ ) , = I t^t - t^t ) 7 ^„ = w - &*%> - ^M a ± Wb7 _1 > = w - Auiu - k x u a , 



(153) 

where / Q are one dimensional distribution functions of the components a. 

This dielectric tensor reduces to the dielectric tensor for plasma in straight magnetic fields for u a = 0. 
It takes a correct account of the Cherenkov-curvature emission and gives the drift corrections to the 
cyclotron emission. We also note that this dielectric tensor is nonhermitian, since k r is not a Killing vector. 
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9. Waves in Cylindrical Coordinates in Anisotropic Dielectric 



In this section we consider polarization of waves in anisotropic dielectric in cylindrical coordinates. 
Since the dielectric properties of a medium are determined by the nonresonant wave-particle interaction, 



we assume that it can be considered in the WKB approximation, so we can use the dielectric tensor (152) 
The properties of a inhomogeneous medium are determined by the very strong, circular magnetic field. 



9.1. Infinitely Strong Magnetic Field 

In the infinitely strong magnetic field the dielectric tensor is (f77j). The dispersion relation is then 



((l-K)(l-nl)-n 2 r ~nl) (1 



n l -nl) = 



(154) 



where — kic/uj. Eq. (154) has solutions 



,(*) 



-{-n x ,0,n r }; 



ni = l- 



1 -K' 



y/(l-K) 2 +Kri{ 



x \ ^/(l-K)(l-K-ni)^-, -n ± , J {I K )(1 - K n\)^ 



(155) 



where e(*) and e"*) are polarization vectors of the t and It mode according to classification of Kazbegi et 



al. 1991. The electric field of the i-wave is always perpendicular to the magnetic field and the wave vector. 
In the limit K — ► the polarization vector for the It wave reduces to 



Ait) 



1 



(156) 



9.2. Finite Magnetic Field 



In the finite magnetic field the dielectric tensor can be found from Eq. ( 152 ). For a cold plasma in the 
center of momentum frame the dielectric tensor is 



1 + 2d 
1 - K 
1 + 2d 



where d = cjp/cj%. Eq. (158). 
Dispersion relation is then 



(-(1 - K)n\ + (1 + 2d)(l - K - nl - n%)) {l + 2d 



n l ~ n l) = o 



where d = ujp/ujg. Eq. ( p.581) has solutions 



n 2 = (1 + 2d), 



e (4) = — {~n x ,0,n r } 
n± 



(157) 



(158) 
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»5 = (i+ag (i-A). 



e a*) 



{ J(l - - X - ^i)^, -n±, — } + 0(d) (159) 
v n i Hi 



v/(l-/f) 2 + ^i V n -L "X 



So that the polarization vectors are the same as in the case of infinitely strong magnetic field within factors 
uj%/ui%. 



9.3. Growth rate of the Cherenkov-drift instability 



Next we calculate the growth rate of the Cherenkov-drift instability ( Lyutikov, Machabeli & Blandford 



1998) of a beam propagating through a medium. In the plane- wave approximation the resonance condition 



is 

1 - k"i 1 / k x Ud \ , 

2 

where we used v res = c(l— ^4-) and e — 1 + 2(5, 5 <C 1. From Eq.( |l6C| ) it is clear, that Cherenkov-drift 

resonance can be satisfied only for superluminal particles with 7 > l/\/25. The emission geometry at the 
Cherenkov-drift resonance is shown in Figs. || and ||. 

Growth rate for the lt-mode in the limit K « is 

u 4ttV f (k^k x u d v<f,k±\ df(p$) 

and growth rate for the t-mode is 

t 4ir 2 q 2 f fk r u d \ 2 df{ P( p) 
= ~m~ J 4 ' V"fcj~c J dp I ^ 4 ^ ~ xUd > ' ^ ' 



It is clear from (161) and (|l62J) that the growth rate of the t-wave is proportional to the drift velocity 
and becomes zero in the limit of vanishing drift velocity. As for the lt-wave, it can be excited in the limit 
of vanishing drift by the conventional Cherenkov mechanism which does not rely on the curvature of the 
magnetic field lines. We recall, that in the limit of a strong magnetic field and oblique propagation lt-wave 



has two branches: one superluminous and one subluminous (Arons & Barnard 1986). On the conventional 



Cherenkov resonance it is possible to excite only subluminous waves. 



10. Conclusion 

In this work we investigated electromagnetic processes associated with a charged particle moving in 
a dielectric in a strong circular magnetic field. We derived a simple expression for the growth rate of 
the Cherenkov-drift instability that, we believe, may be responsible for the generation of the pulsar radio 
emission. We leave the discussion of the astrophysical applications of our results for the subsequent paper 
(Lyutikov, Blandford & Machabeli 1998). Here we just note, that the Cherenkov-drift developing on the 
open field lines in the outer parts of pulsar magnetosphere can explain various features of the cone-type 
(Rankin 1992) emission patterns observed in pulsars. 
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Fig. 1. — Geometry of the considered problem. The magnetic field lines are concentric coplanar circular 
arcs with the radius of curvature much larger than the size of the region. 
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Fig. 2. — Surfaces of constant phase of Hankcl function H u ' for v = 5. Circle Xr = v is an analog of a light 
cylinder. For radii much smaller than the "light cylinder" radius r v = v/X, the rotation is similar to solid 
body rotation, while for the radii much larger than the light cylinder radius the surface of a constant phase 
has a form of unwinding spiral with a wavelength Xrjv. 
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Fig. 3. — Polarization of normal modes in the limit of very strong magnetic field. The electric field vector 
of the t-mode is in the plane e r — e x and the electric vector of the ^-mode is orthogonal to e t — k plane. 
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Fig. 4. — Cherenkov-drift emission in the case S <C u 2 d /c 2 . Drift velocity Ud is perpendicular to the plane of 
the curved field line (B — R c plane, R c is a local radius of curvature). The emitted electromagnetic waves 
are polarized along u^. The emission is generated in the cone centered at the angle 9 em = Ud/c and with 
the opening angle (2(5) 1 I 2 < 6 em . 
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Fig. 5.— 



Emission geometry of the Chcrcnkov-drift resonance 



